Non-Gaussian quantum discord for Gaussian states 
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In recent years the paradigm based on entanglement as the unique measure of quantum correla- 
tions has been challenged by the rise of new correlation concepts, such as quantum discord, able to 
reveal quantum correlations that are present in separable states. It is in general difficult to compute 
quantum discord, because it involves a minimization over all possible local measurements in a bi- 
partition. In the realm of continuous variable (CV) systems, a Gaussian version of quantum discord 
has been put forward upon restricting to Gaussian measurements. It is natural to ask whether 
non- Gaussian measurements can lead to a stronger minimization than Gaussian ones. Here we fo- 
cus on two relevant classes of two-mode Gaussian states: squeezed thermal states (STS) and mixed 
thermal states (MTS), and allow for a range of experimentally feasible non-Gaussian measurements, 
comparing the results with the case of Gaussian measurements. We provide evidence that Gaussian 
measurements are optimal for Gaussian states. 
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I. INTRODUCTION 

In recent years the paradigm based on entanglement [l[ 
as the unique genuine measure of quantum correlations 
has been challenged by the argument that the notion 
of nonseparability may be insufficient to encompass all 
correlations that can be fairly regarded as quantum, or 
nonclassical. This has given spur to the development 
of conceptually new correlation measures, such as quan- 
tum discord 0tJ1 > based on local measurements and able 
to reveal quantum correlations that are present even in 
separable states. These correlations can be interpreted 
as an extra amount of information that only coherent 
operations can unlock [5j]. In fact, there are several in- 
dications suggesting that general quantum correlations 
might be exploited in quantum protocols @, including 
mixed state quantum computation Q and remote state 
preparation [§]. Therefore, a more complete theoretical 
and experimental investigation thereof is now a central 
issue in quantum science and technology @. 

The definition of discord involves an optimization over 
all possible local measurements in a bipartion, the opti- 
mal measurement leading to a minimal value of quantum 
discord. To perform the optimization is remarkably diffi- 
cult, which hampers analytical progress in the area. This 
fact has led to the definition of other correlation measures 
which are conceptually similar but easier to compute, 
such as the geometric discord [lOj . In the realm of finite- 
dimensional systems, where the concept of discord was 
first introduced, analytic results for quantum (geomet- 
ric) discord have been obtained for pairs of qubits when 
the global state is in X form (in arbitrary form) [l(| [H| ■ 

In the realm of continuous variable (CV) systems, ini- 
tial research efforts on quantum discord have focused on 
Gaussian measurements. The Gaussian quantum dis- 
cord, proposed in [l2|, HH, is defined by restricting the 



minimization involved in the definition of discord to the 
set of Gaussian POVMs [14j and it can be analytically 
computed for Gaussian states. Its behavior in noisy chan- 
nels has been studied in Ref. [l5| - where it was shown 
that it is more robust than entanglement to the decor- 
relating effect of independent baths and more likely to 
yield non-zero asympotic values in the case of a common 
bath - while its relation to the synchronization proper- 
ties of detuned, correlated oscillators has been analysed 
in Ref. [3. 

It is natural to investigate CV quantum discord beyond 
Gaussian measurements: non-Gaussian ones may indeed 
allow for a stronger minimization of discord, and in this 
case the Gaussian discord would be an overestimation of 
the true discord. Here we focus on Gaussian states and 
ask whether Gaussian measurements are optimal in this 
case, i.e., whether the Gaussian discord is the true dis- 
cord for Gaussian states. This question is relevant for 
two main reasons: On one hand, if discord is a truly 
useful resource for quantum information protocols @,|g], 
then it is crucial to have a reliable estimate of its actual 
value. On the other hand, from a fundamental point of 
view it is important to establish how different kinds of 
measurements can affect correlations in quantum states. 
A further motivation comes from the fact that indeed 
for some non-Gaussian states e.g., CV Werner states, 
non-Gaussian measurements such as photon counting has 
been proven to lead to a better minimization jl7l |. 

The optimality of Gaussian measurements has already 
been proven analytically for two-mode Gaussian states 
having one vacuum normal mode (13(, by use of the so- 
called Koashi- Winter relation [l8[ , but no analytic argu- 
ment is available in the general case. We address the 
question numerically, for the case of two-modes, upon 
considering two large classes of Gaussian states, the 
squeezed thermal states (STS) and the mixed thermal 
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states (MTS), and allowing for a range of experimentally 
feasible non-Gaussian measurements based on orthogo- 
nal bases: the number basis, the squeezed number basis, 
the displaced number basis. As a result, we provide evi- 
dence that Gaussian quantum discord is indeed optimal 
for the states under study. In addition, we also investi- 
gate the CV geometric discord [19|, comparing the case 
of Gaussian and non-Gaussian measurements. 

This work is structured as follows. In sec. [IT] we review 
quantum discord and the Gaussian version of it; in sec. 
IIIII we thoroughly describe the basic question we want to 
address in this work and introduce non-Gaussian mea- 
surements and non-Gaussian discord; in sec. QVl E] I VII 
we present our key results concerning non-Gaussian dis- 
cord upon measurements in the number basis, squeezed 
number basis and displaced number basis; in sec lVIII we 
discuss the behavior of non-Gaussian geometric discord; 
finally, sec. I Villi closes the paper discussing our main 
conclusions. 



II. QUANTUM DISCORD AND GAUSSIAN 
DISCORD 

Starting from the seminal works by Ollivier and 
Zurek [2fl and Henderson and Vedral [3], various mea- 
sures of quantum correlations which go beyond the tra- 
ditional entanglement picture have been defined Q . The 
most common measure of such correlations is the quan- 
tum discord 0, H[. Let us consider a bipartite system 
composed of subsystems A and B. The total correlations 
in the global state are measured by the mutual informa- 
tion I(A : B) = S{q a ) + S{q b ) - S{g AB ). Whenever 
I {A : B) > 0, the subsystems are correlated and we can 
gain some information about A by measurements on B 
only. However, there is no unique way of locally probing 
the state of B: to do it, we can perform different local 
measurements, or POVMs. Any such local POVM is 
specified by a set of positive operators {H B = M^Mp^} 
on subsystem B summing up to the identity, = ^- 

When measurement result x is obtained, the state of A is 
projected onto g x A = Tyb[M b q A bM' b ]. The uncertainty 
on the state of A before the measurement on B is given 
by S(qa), while the average uncertainty on the state of A 
after the measurement is given by the average conditional 
entropy S n (A\B) = J2 x PxS(g x A ). Their difference 

S(q a ) - S n (A\B) = S(g A ) - ]T^S(^) 

X 

represents the average gain of information about the state 
of A acquired through a local measurement on B. The 
maximal gain of information that can be obtained with 
a POVM, 

C(A : B) = max \S(g A ) - S n (A\B)} = (1) 

{YldPOVM} 

= S{q a ) - min{nePOVM}[S U (A\B)} 



coincides with the measure of classical correlations orig- 
inally derived in Q under some basic and natural re- 
quirements for such a measure. Quantum discord is then 
defined as the difference between the mutual information 
and the classical correlations: 

D(A : B) = I(A : B) - C(A : B) (2) 

and measures the part of the total correlations that can- 
not be exploited to gain information on A by a local 
measurement on B, i.e., measures the additional quan- 
tum correlations beyond the classical ones. 

It can be verified (see e.g. [l(|) that the classical corre- 
lations coincide with the mutual information in the sys- 
tem after the measurement, maximized over all possible 
POVMs: 

C(A:B)=ma HnePOVM} I n (A:B) (3) 

where I n {A : B) = S(q%) + S(gg) - S{g%) and 
the unconditional post-measurement states are given by 
Q% = E X M BQAB M X J, & = Tr s [£ x M%g AB M B % 
= Tr A [Z, x M B QabM^}. Therefore, the quantum 
discord coincides with the difference between the mu- 
tual information before and after the measurement, min- 
imized over all possible POVMs: 

D(A : B) = mm {nePOVM} [I(A : B) - I n (A : B)] (4) 

From the prevoius considerations, it is clear that D(A : 
B) = if and only if there is a local measurement n# 
which leaves the global state of the system unaffected: 
3n, q A b — Q^b- Such states are called quantum- classical 
states and are in the form 

Xab = ^PiQA,i ® \i){A (5) 

i 

where pi is a probability distribution and {\i)} is a basis 
for the Hilbert space of subsystem B. For such states, 
there exists at least one local measurement that leaves 
the state invariant and we have I (A : B) = C(A : B), 
which means that we can obtain maximal information 
about subsystem A by a local measurement on B without 
altering the correlations with the rest of the system. 

In the realm of continuous- variable systems, the Gaus- 
sian discord [l2l [l3j is defined by restricting the set of 
possible measurements in Eq. (UJ to the set of Gaussian 
POVMs and minimizing only over this set. The 

Gaussian discord can be analytically evaluated for two- 
mode Gaussian states, where one mode is probed through 
(single-mode) Gaussian POVMs. The latter can be writ- 
ten in general as 

n B (ri) =ir- 1 D B {ri)Q M tf B {n) 

where Db(t]) = exp(?7&^ — 77*6) is the displacement oper- 
ator, and qm is a single-mode Gaussian state with zero 



mean and covariance matrix om 




. Two-mode 
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Gaussian states can be characterized by their covariance 



matrix 



ix uab — ^ qt ^ ^ ■ By means of local 



unitaries 



that preserve the Gaussian character of the state, i.e. lo- 
cal symplectic operations, o ab may be brought to the so- 
called standard form, i.e. A — diag(a, a), B — diag(6, b), 
C = diag(ci,C2). The quantities I\ — detA, I2 — detB, 
I3 = dot C, I4 = det cab are left unchanged by the trans- 
formations, and are thus referred to as symplectic invari- 
ants. The local invariance of the discord has therefore 
two main consequences. On the one hand, correlation 
measures may be written in terms of symplectic invari- 
ants only. On the other hand, we can restrict to states 
with o already in the standard form. Before the mea- 
surement we have 

S{ QAB ) = h(d+) + h(d-), (6) 
S(qa) = h{y/h), S(g B ) = HVh) (7) 

where h[x] = (x+1/2) log(x+l/2)-(x-l/2) log(a;-l/2) 
and d± are the symplectic eigenvalues of qab expressed 
by d\ = 1/2[A± V / A^ -4/4, A = h+I 2 + 2h. After the 
measurement, the (conditional) post- measurement state 
of mode A is a Gaussian state with covariance matrix 
op that is independent of the measurement outcome and 
is given by the Schur complement op = A — C(B + 
o'm)~ 1 C t . The Gaussian discord is therefore expressed 
by 



D y (A : B) = hWL 



h(d-)-h(d+) 
V M /i(det y/op) 



(8) 



where we use two key properties: i) the entropy of a 
Gaussian state depends only on the covariance matrix, 
and ii) the covariance matrix op of the conditional state 
does not depend on the outcome of the measurement. 
The minimization over om can be done analytically. For 
the relevant case of states with C — diag(c, ±c), includ- 
ing STS and MTS (see below), the minimum is obtained 
for a — f3 — 1/2, 7 = i.e. when the covariance matrix of 
the measurement is the identity. This corresponds to the 
coherent state POVM, i.e. to the joint measurement of 
canonical operators, say position and momentum, which 
may be realized on the radiation field by means of hetero- 
dyne detection. For separable states the Gaussian discord 
grows with the total energy of the state and it is bounded, 
D < 1; furthermore, we have D — iff the Gaussian state 
is in product form qab — QA® Qb- 



III. NON GAUSSIAN DISCORD 

In this work we consider Gaussian states, and ask 
whether non-Gaussian measurements can allow for a bet- 
ter extraction of information than Gaussian ones, hence 
leading to lower values of discord. 

The optimality of Gaussian measurements has been 
already proven for a special case 13]: that of two-mode 



Gaussian states having one vacuum normal mode. Indeed 
any bipartite state qab can be purified, qab ==> \"4>)abc'-> 
then, the Koashi- Winter relation 1181. 



D(A : B) = Ej(A : C) + S(g B ) ~ S(q A b) 



(9) 



relates the quantum discord D and the entanglement of 
formation Ef of reduced states qab and Qac respectively. 
Given a (mixed) two-mode Gaussian state qab, there ex- 
ists a Gaussian purification \ip) abc- In general, the pu- 
rification of qab requires two additional modes, so that 
qac is a three-mode Gaussian state. In the special case 
when one normal mode is the vacuum, the purification 
requires one mode only. In this case, qac represents a 
two-mode Gaussian state and Ef(A : C) can be evalu- 
ated [H. From E f (A : C), by means of Eq. ([3), one 
can obtain D(A : B) (the exact discord) and a compari- 
son with D G (A : B) proves that D(A : B) = D g (A : B). 

In the general case, there is no straightforward ana- 
lytical way to prove that Gaussian discord is optimal. 
Therefore, we perform a numerical study. Since tak- 
ing into account the most general set of non-Gaussian 
measurements is an extremely challenging task, one can 
rather focus on a restricted subset. We choose to focus 
on a class of measurements that are realizable with cur- 
rent or foreseable quantum optical technology. These are 
the the projective POVMs, II = {IT n }, represented by 
the following orthogonal measurement bases: 



H„ = D(a)S(r)\n)(n\S( r y D( a y , n = 



,00 (10) 



and D(a) 



where S(r) — exp (— r*^- , „ 
exp(aa^ — a* a) are respectively the single- mode squeez- 
ing and displacement operators [2lJ. The set of projec- 
tors in (fTUf is a POVM for any fixed value of a and r. If 
a = r = we have the spectral measure of the number 
operator, describing ideal photon counting IT n = |n)(n|. 
If a > 0,r = we are projecting onto displaced num- 
ber states l2p . if a = 0, r > onto squeezed number 
states [23l426| . While more general non Gaussian mea- 
surements are in principle possible, the class (|10[) encom- 
passes most of the measurements that can be realistically 
accessed. 

In the following, we will evaluate the non-Gaussian 
quantum discord defined by 



D Aiy (A : B) = hWh) - h(d-) - h{d+) + S n '* y (A\B) 

(11) 

where the non- Gaussian measurements are given by Eq. 
(fTUf above. For the non-Gaussian conditional entropy we 
have 

5 n ^ e ( J 4| J B)=^p„5(^,„), 



QA.n = — Tr s [n n gAsII n ] 

Pn 

p n = Tr AB[HnQABHn] 



(12) 
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In the following we consider two classes of Gaussian states 
in order to assess the performances of the above measure- 
ments. These are the two-mode squeezed thermal states 
(STS) H2HI3: 

g = S(\)iy 1 (N 1 )®^ 1 (N 2 )S(\^ (13) 

and the two-mode mixed thermal states (MTS) [3p| 

p = U{4>)v 1 {N 1 )®v 1 {N 2 )U{^ (14) 

where i^(iVj) are 1-mode thermal states with thermal 
photon number Ni) S(X) = exp{A(a|a2 — 0102)} is the 
two-mode squeezing operator (usually realized on optical 
modes through parametric down-conversion in a nonlin- 
ear crystal); and U{<f>) — exp{(f>(a\a 2 — a\a 2 )} is the two- 
mode mixing operator (usually realized on optical modes 
through a beam splitter). 

In particular, in the following we will focus on the sim- 
plest case of symmetric STS with N\ = N 2 G [10~ 5 ,1] 
A G [0, 0.5]. As for MTS, we cannot consider the symmet- 
ric case (since if N% = N 2 then the mutual information 
vanishes and there are no correlations in the system), 
therefore we consider the unbalanced case and focus on 
G [0,7r/2] and Nx,N 2 e [10~ 5 ,1]. 

IV. NUMBER BASIS 

Let n„ = |n)(n|. In this case, the post-measurement 
state is 

Qn ® \ n )( n \ = S(h,k),(n,n)\h){k\ ® \n){n\ (15) 

\h,k J 

and we have the following expression for the density ma- 
trix elements 

Q (hlkUn , n) =Y,PT( N M h (N2)O hn (st)0* kn (st) (16) 

s,t 

where pf(N) = N s (1 + N)'^ s+1 ^ and O hn (st) = 
(hn\0\st) with O = S(\),U(<j>) for STS and MTS re- 
spectively. The post-measurement state g^ is diagonal 
(see appendix E|, 

( h \6n\ k ) = S hk Q(h,h),(n,n) (17) 

As a consequence, the entropy of the post-measurement 
state can be expressed as: S(g^) — H({g^ hh ^ ( n ,n)}) — 
H(p(A\B = n)) where H is the Shannon en- 
tropy of the conditional probability p(A\B = n) = 
(p(0, n),p(l, n), ■ ■ ■ )/p n , and therefore the overall condi- 
tional entropy can be simply expressed in terms of the 
photon number statistics: 

S(A\{U n }) = Y J PnH{p{A\B = n)) = 

n 

= H(p(A,B))~H(p(B)) (18) 




FIG. 1: Gaussian and non-Gaussian quantum discord for 
STS as a function of A, for different values of Ni = N2 



with p(A, B) = {p{A = n,B = m)} and p(B) = {p(B = 
n)}. In view of this relation, the only elements of the 
number basis representation of the density matrix g that 
are needed are the diagonal ones, i.e. one has to de- 
termine the photon number statistics for the two-mode 
STS or MTS state. The required matrix elements can 
be obtained in terms of the elements of the two-mode 
squeezing and mixing operators (see appendix |A| . One 
has of course to define a cutoff on the dimension of the 
density matrix. This can be done upon requiring that the 
error on the trace of each state considered be sufficiently 
small: e err = 1 — Trg < 10 -3 . 

We have compared Gaussian and non-Gaussian quan- 
tum discord (with the non-Gaussian measurements corre- 
sponding to photon number measurements) for STS and 
MTS states with a wide range of squeezing, mixing and 
thermal parameters. In Fig. [T] we show results for STS 
with varying A and N x = JV 2 = 10~ 2 , Nx = N 2 = 1. The 
key result is that the non-Gaussian quantum discord is 
always greater than its Gaussian counterpart for all val- 
ues of Ni and A. The gap grows with increasing N± and 
A. In Fig. |we show results for MTS Ni = {0.1, 1} and 
q = N 2 /Ni = {0, 0.1, 0.4, 0.5}. Also in this case, the non- 
Gaussian discord is always higher than the Gaussian one. 
Both results indicate that the Gaussian (heterodyne) 
measurement is optimal for STS and MTS states, at least 
compared to photon counting, in the sense that it allows 
for a better extraction of information on mode A by a 
measurement on mode B. 



V. SQUEEZED NUMBER BASIS 

We now analyze the case of non- Gaussian mea- 
surements represented by the squeezed number ba- 
sis |?v)(n r | = S(r)\n){n\S(r)\ where S(r) = 

exp (— r*^ — r ^ a 2 ) is the single mode squeezing oper- 
ator. A local measurement in the squeezed number ba- 
sis is equivalent to a measurement in the number basis, 
performed on a locally squeezed state. In formulas, the 
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FIG. 2: Gaussian and non-Gaussian quantum discord for 
MTS states as a function of 4> f° r different values of Ni and 

q = JVa/JVi 




FIG. 3: Gaussian and non-Gaussian quantum discord for STS 
with Ni = 1 as a function of A and for different values of local 
squeezing r 



probability of measuring n r on one subsystem when the 
state is the g is 

Pg (n r ) = Tr(l <g> \n r )(n r \g) = Tr(l <g> \n)(n\S\r)gS(r)) = 
= Tr(l® |n)(n|g r ) = p er (n) (19) 

i.e., is equal to the probability of measuring n on the 
locally squeezed state g r = S(r)gS{r)\ and the relative 
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FIG. 4: Gaussian and non-Gaussian quantum discord for 
MTS states for Ni = 1, N2 = as a function of 4> an d for 
different values of local squeezing r 



post-measurement state is 

gf lr = Tr B [H ® \n r )(n r \gt ® \n r ){n r \]/p e (n r ) = 

= Tr s [l ® \n)(n\g r t 8 \n)(n\]/ P r g (n) = gf n (20) 

The general idea is that measurements on a state g in a 
basis that is obtained by performing a unitary (Gaus- 
sian) operation V on the number basis \n)(n\ can be 
represented as measurements on the number basis of a 
modified state gv — V gV^ on which the local unitary 
operation acts. 

In the case of the squeezed number basis, the post- 
measurement state is not diagonal, therefore the rea- 
soning leading to Eq. (|18l) does not hold. The 
post-measurement state matrix elements (g^ n )h,k = 
Q{h.k),{n,n) can be obtained directly by evaluating the 
expression (|16[) where now the expression Ohk{st) — 
(hk\0\st) (where O = S(X),U((f>)) must be substituted 
with 0' hk (st) = (hk\S(r)0\st) = E q (k\S(r)\q)(hq\0\st), 
and the elements of the single mode squeezing operator 
are given in [3l[ (eq. 20) or in [32[ (eq. 5.1). 
We have evaluated the Gaussian and non-Gaussian quan- 
tum discord for STS and MTS states with a wide range 
of two-mode squeezing and thermal parameters. Non- 
Gaussian measurements are done in the squeezed pho- 
ton number basis, II n = S(r)\n)(n\S(rY with variable 
r G [0,0.5]. The effect of local squeezing on non- 
Gaussian quantum discord is negligible in the whole pa- 
rameter range under consideration: we compare the non- 
Gaussian discord for different values of r and find that 
all curves collapse. This can be seen in fig. [3] and fig. 
H where plot the behavior for N ± = N 2 = 0.01 (STS) 
and Ni = 1,N 2 = (MTS). The same behavior is ob- 
served in the whole parameter range under investigation. 
We have verified numerically that the post-measurement 
states of mode A gfi are not equal as r varies (i.e., the 
post-measurement states corresponding to measurement 
result n r change with r), yet the sum ^2 n p n S(g^ n ) is 
equal for all values of r under investigation. Therefore, 
the squeezing in the measurement basis has no effect on 
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FIG. 5: Gaussian and non-Gaussian quantum discord for STS 
with Ni = 1 as a function of A and for different values of local 
displacement a 



FIG. 6: Gaussian and non-Gaussian quantum discord for 
MTS states for Ni = 1, N2 = as a function of 4> an d for 
different values of local displacement a 



the discord, at least for the values of squeezing consid- 
ered: in particular, it cannot afford a deeper minimiza- 
tion than that obtained without local squeezing. This 
indicates that the heterodyne measurement remains op- 
timal also with respect to measurement in the squeezed 
number basis. 



VI. DISPLACED NUMBER BASIS 

We finally analyze the case of non-Gaussian mea- 
surements represented by the displaced number basis 



exp(aa 



\n a )(n a \ = D(a)\n){n\D{a)\ where D(a) 
a* a) is the single mode displacement operator. Accord- 
ing to the general considerations above, a local measure- 
ment in the displaced number basis is equivalent to a 
measurement in the number basis, performed on a locally 
displaced state g a . As in the case of the squeezed number 
basis, the post-measurement state is not diagonal and we 
need all matrix elements (Qa n )h,k = Q(h,k),(n,n)- They 
can be obtained directly by evaluating the expression 
(fT6|) where the expression Ohk(st) = (hk\S(X)\st) (where 
O = S(X),U(4>)) must be substituted with 0' hk \st) = 
(hk\D(a)0\st) = ^ q (k\D(a)\q)(hq\0\st), and the ele- 
ments of the single mode displacement operator are given 
in [13 (eq. 1.46). 

The evaluation of the non-Gaussian quantum discord 
can be simplified by first noticing that one can consider 
real values of a only. Indeed, the quantum discord only 
depends on the modulus \a\. This is shown in detail 
in the appendix [B] by using the characteristic function 
formalism. Consider g^ , the post-measurement state of 
mode A after measurement result n a is obtained on B. 
If we change the phase of a, a —> a' = e l6 a we find that 



(21) 



where U is a unitary operation corresponding to a simple 
quadrature rotation 



«i 



die 



,1) 



Therefore, we have g^, 7^ g£ , but g^, and gfc have 
the same spectrum, since they are related by a unitary. 
Therefore, the entropy of the reduced post-measurement 
state g£ does not depend on the phase of a but just on 
\a\. If follows that the non-Gaussian quantum discord of 
g a does not depend on the phase of a. 
We have evaluated the Gaussian and non-Gaussian quan- 
tum discord for STS and MTS states with a wide range 
of two- mode squeezing and thermal parameters. Non- 
Gaussian measurements are done in the displaced pho- 
ton number basis, II n = D(a)\n)(n\D(ay with variable 
a £ [0,2.5]. In fig. [5] and fig. [5] we plot the Gaussian 
and non-Gaussian quantum discord. We see that greater 
displacements lead to lower values of the non-Gaussian 
quantum discord, but the decrease is insufficient to match 
the Gaussian quantum discord, which remains optimal. 
However, the non-Gaussian quantum discord approxi- 
mates the Gaussian one as a — s- 00. This is analytically 
proven below in the appendix [Cl There we find that for 
both STS and MTS 



Oil 



as a 



00 



(23) 



at e" 



,0 



(22) 



i.e, the conditional states g£ becomes independent of 
n and equal to the n = result. As a consequence, 
the conditional entropy in the displaced number basis is 
equal to the entropy of the post-measurement state for 
any measurement result, and, in particular, for n = 0: 

S n ^(A\B)=Y,P" S (et)^S(ei) as a ^ co 

n 

(24) 

But g^ Q is just the post-measurement state we obtain 
after a heterodyne detection on mode B (equal for all 
measurement result modulo a phase space translation 
which is irrelevant as for the entropy). Therefore, we 
also have S°(A\B) = S(g^ ) and the non Gaussian dis- 
cord D-^^ (A : B) in the displaced number basis tends to 
the Gaussian discord D s (A : B) as a — !• 00. 

Actually, we cannot prove that the D^ g (A : B) is 
lower bounded by D s (A : B), and we cannot rule out 
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the possibility that D NS {A : B) < D G {A : B) for in- 
termediate values of a. However, our numerical data 
do not support this possibility since we never observe 
D MG {A : B) < D G (A: B) and we expect that D N9 {A : 
B) — > D G (A : B) from above as a — > oo. 
In conclusion, we have analytical and numerical evidence 
that the heterodyne measurement remains optimal also 
with respect to measurement in the displaced number 
basis. 



VII. GEOMETRIC DISCORD 

In this section, we briefly consider the recently intro- 
duced measure of geometric discord and compare results 
with those obtained for the quantum discord. The geo- 
metric discord is defined as 

D G (A : B) = min \\g AB - X ab\\2 (25) 
{xA B ec} 

and it measures the distance of a state from the set C 
of quantum-classical states where \ \A\\2 = Tr[A^ A] is the 
Hilbert-Schmidt distance. Clearly Dq — iff D = 0, 
since both measures vanish on the set of classically cor- 
related states. In particular, it has been be proven that 
Dq can be seen a measure of the discrepancy between a 
state before and after a local measurement on subsystem 
B 0: 

D G (A:B)= min \\q ab - Qab\\ 2 (26) 
{tigpovm} 

where the unconditional post-measurement state is given 
by g% = J2x M b8abM^ . Notice that D G and D are 
not monotonic functions of one another and the relation 
between them is still an open question. However, in many 
cases Dq is much simpler to evaluate than D. 

Analogous to the case of Gaussian discord, a Gaussian 
version of the geometric discord can be defined by re- 
stricting to Gaussian measurements [19]. Again, it can 
be analytically computed for two-mode Gaussian states. 
With the same reasoning of sec. |TT]one easily obtains 

D%{qab) = min . M Tr[(g J 4B - Qp ® Qm) 2 ] (27) 

Exploiting the property that Tr[gi£2] = l/dct[(o"i + 
<72)/2], for any two Gaussian states g\ and Qi, 

Dq (A: B) — min CTjv/ {1/Vdet cab 4(28) 
+l/v / det(cr P © om) - 2/Vdet [(ctab + a P © er M )/2]} 

For for the relevant case of STS and MTS, the minimum 
is obtained with the <jm elements given by a — f3 = 

v/ ^ (V4afc 3 - a 3c2+v ^ , 7 = 0. The least disturbing Gaus- 
sian POVM for STS, according to the Hilbert-Schmidt 
distance, is thus a (noisy) heterodyne detection, a result 
which is analogous to what found in the case of quantum 
discord. If one constrains the mean energy per mode, the 
Gaussian quantum discord gives upper and lower bounds 
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FIG. 7: (Top) Gaussian and non-Gaussian geometric discord 
for SFS with iVi = 1 as a function of A and for different values 
of local squeezing r; (Bottom) Gaussian and non-Gaussian 
geometric discord for MTS states for Ni = 1, JV2 = as a 
function of cj> and for different values of local squeezing r 



to the Gaussian geometric discord. In absence of such a 
provision, the geometric discord can vanish for arbitrar- 
ily strongly nonclassical (entangled) Gaussian states, as 
a consequence of the geometry of CV state spaces. 

Also in this case, we may consider non-Gaussian mea- 
surements and evaluate a non-Gaussian geometric dis- 
cord: 

D^(A:B)=Tr[(g AB ~g%) 2 } (29) 

For measurement in the number basis, we can easily ob- 
tain 

D^ =n( Q )+J2\(pn\Q\qn)\ 2 (30) 

npq 

where ll(q) = — j= is the (Gaussian) state pu- 

4-\/det(er) 

rity [33| . In the case of measurements in the squeezed or 
displaced number basis, we have to use g r and g a instead 
of q in Eq. (l30l) . In general, in order to compute the 
geometric discord we need to compute matrix elements, 
and we use the same numerical methods described above. 
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FIG. 8: (Top) Gaussian and non-Gaussian geometric discord 
for STS with Ni — 1 as a function of A and for different values 
of local displacement a (Bottom) Gaussian and non- Gaussian 
geometric discord for MTS states for Ni = 1, N2 = as a 
function of <j> and for different values of local displacement a 



A. Results 

. We have compared the Gaussian and non-Gaussian 
geometric discord for STS and MTS in a wide range of 
parameters. We have considered measurements in the 
number, squeezed number and displaced number basis 
for the same values of the parameters given in the pre- 
ceding sections. Results are plotted in Figs. [7] and [5] In 
general, at variance with the results for quantum discord, 
we find that non-Gaussian measurements can provide 
lower values of geometric discord than Gaussian ones. 
Among the class of non-Gaussian measurements we have 
considered, the optimal one is provided by the number 
basis, which gives values of geometric discord that are 
always lower than those given by the optimal Gaussian 
measurement. The non-Gaussian geometric discord in- 
creases for increasing r and a, and it can become greater 
than its Gaussian counterpart. These results are very 
different from the quantum discord case: on one hand, 
the (non-Gaussian) geometric discord is substantially af- 
fected by the local squeezing; on the other hand, it does 
not approach the Gaussian one when the displacement 
a — > 00, but it grows monotonically. Indeed if we increase 
the squeezing or displacement in the measurement basis, 
the post-measurement state is more distant (in Hilbert- 



Schmidt norm) from the original one. As already noticed, 
performing the measurement is the squeezed (displaced) 
number basis in equivalent to first squeezing (displacing) 
the state and then measuring it in the number basis. 
The local squeezing and displacement have the effect of 
increasing the energy of the state, shifting the photon 
number distribution P(B = n) towards greater values 
of n. This causes the overlap between the post mea- 
surement state and the original state to decrease, and 
therefore their distance to increase. 

Let us futher comment on the difference between the 
quantum discord and the geometric discord cases. Quan- 
tum discord and geometric discord both vanish for clas- 
sical states, but are not monotonic functions of one an- 
other, and thus they are truly different quantities. The 
geometric discord, based on the Hilbert-Schmidt dis- 
tance, is a geometric measure of how much a state is per- 
turbed by a local measurement, while quantum discord 
assesses to which extent correlations are affected by a 
local measurement. While for the quantum discord well- 
defined operational and informational interpretations can 
be found d, @ , for the geometric discord the situation is 
more problematic. Indeed, one can design protocols in 
which the geometric discord can in some cases be related 
to the protocols' performances @, |35|; however, recent 
discussions [36| . show that, as consequence of the non- 
invariance of the Hilbert-Schmidt norm under quantum 
evolutions, it is difficult to find a conclusive argument 
about the relevance of geometric discord as a measure of 
quantumness of correlations. Our data show that non- 
Gaussian measurements can yield optimal values of the 
geometric discord, contrary to the case of quantum dis- 
cord. Hence, the behavior of quantum discord and ge- 
ometric discord with respect to different types of mea- 
surements is different. This is a further indication that 
the geometric discord cannot be used as a good bench- 
mark for the quantum discord and that the degree of 
quantumness measured, if any, by such a quantity has a 
fundamentally different nature. 



VIII. DISCUSSION AND CONCLUSIONS 

The definition of discord involves an optimization over 
all possible local measurements (POVMs) on one of the 
subsystems of a bipartite composite quantum system. In 
the realm of continuous variables (CV), initial research 
efforts on quantum discord restricted the minimization 
to the set of (one- mode) Gaussian measurements. 

In this work we have investigated CV quantum dis- 
cord beyond this restriction. We have focused on Gaus- 
sian states, asking whether Gaussian measurements are 
optimal in this case, i.e., whether the Gaussian discord 
is the true discord for Gaussian states. While a posi- 
tive answer to this question had already been given for 
the special case of two-mode Gaussian states having one 
vacuum normal mode (by means of an analytical argu- 
ment based on the Koashi- Winter formula), no general 
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result was available so far. We have addressed our central 
question upon considering two large classes of two-mode 
Gaussian states, the squeezed thermal states (STS) and 
the mixed thermal states (MTS), and allowing for a wide 
range of experimentally feasible non-Gaussian measure- 
ments based on orthogonal bases: the photon number 
basis, the squeezed number basis, the displaced number 
basis. For both STS and MTS states, in the range of pa- 
rameters considered, the Gaussian measurements always 
provide optimal values of discord compared to the non- 
Gaussian measurements under analysis. Local squeezing 
of the measurement basis has no appreciable effect on 
correlations, while local displacement leads to lower val- 
ues of the non-Gaussian discord, which approaches the 
Gaussian one in the limit of infinite displacement. 

Overall, for the explored range of states and measure- 
ments, we have evidence that the Gaussian discord is the 
ultimate quantum discord for Gaussian states. We note 
that the optimality of Gaussian measurements suggested 
by our analysis is a property which holds only for Gaus- 
sian states. In the case of non-Gaussian states, e.g., CV 
Werner states, non-Gaussian measurements such as pho- 
ton counting can lead to a better minimization, as was 
recently proven in Ref. [17] . 

Wc also have investigated the CV geometric dis- 
cord [l!|, comparing the Gaussian and non-Gaussian 
cases. We have shown that the behavior of geometric 
discord is completely different from that of quantum dis- 
cord. On one hand, non-Gaussian measurements can lead 
to lower values of the geometric discord, the number basis 
measurement being the optimal one; on the other hand, 
the effects of both local squeezing and displacement are 
strong and consist in a noteworthy increase in the non- 
Gaussian geometric discord. The remarkable differences 
between quantum and geometric discord imply that the 
latter cannot be used as a benchmark of the former. 

Both in the case of the discord and geometric discord a 
definite answer on the optimal measurement minimizing 
the respective formulas would require the extension of 
the set of non-Gaussian measurements to possibly more 
exotic ones and the application of those realizable in ac- 
tual experiments to a broader class of Gaussian and non- 
Gaussian states. While we leave this task for future re- 
search, our results on discord support the conjecture that 
Gaussian measurements are optimal for Gaussian states 
and allow to set, for the class of states analyzed, a tighter 
upper bound on the entanglement of formation for 1x2 
modes Gaussian states, via the Koashi- Winter relation. 



indeed: 

g = Y / P t s h (Ni)p t t h (N 2 )0\st)(st\0^ = 

s,t 

= e ihm [Y.pfpi h °Ust)oi m { S t)Xk2) 

(h,n),(k,m) V s,t J 

where pf{N) = N s (1 + 7V)-( S + 1 ) where O hn (st) = 
(hn\0\st) and 0* krn (st) = (st\0^\km) = (km\0\st)* , 
where O = S(X), O = U(<f>) for STS and MTS respec- 
tively. The post measurement states can be written as: 

ei <g> \n)(n\ = 0(h : k),(n,n)\h)(k\ ® \n)(n\ (A3) 

\h,k J 

and therefore we need to evaluate the matrix elements 
0(M),(n,n) = 52pl h pl h O hn (st)0* kn (st) (A4) 

s,t 

The elements of the two-mode squeezing operator are 
given in [26] (cq. 22): 

(hn\S(X)\st) = 5 t+ h, s +nf X (h,n,s,t) = S t+hiS+n x 

min(s,i) min(h,n) 



E E (-irWA) 



t-\-h—a—b—l a—b+h — s 

fj, X 



a=0 6=0 



{t + h - a - b)\[s\t\h\n^/ 2 
a\(t - a)\(s - a)\b\(n - b)\(h - b)\ 



(A5) 



where /j, = e x , while the elements of the two- mode mixing 
operator 

min{/i,s} 

(hn\U(</>)\st) = 5 h+n , s+t E A ah-a 

a-max{0,/i-i} 

h\{s + t-hy.. h _ a 



— 3h+n,s+t 

a-max{0,/i-t} 

'S\ ( t 



slti 



.(_!)#»-«• x 



a J \ h — a 



sm<b s+h - 2a cos6 t+2a - h 



(A6) 



In order to evaluate (|A4j) . we need Ohn(st)Ol n (st). Due 
to the S's appearing in both (|A5j) and (|A6j) . the following 
relations must be satisfied: 



Appendix A: The post-measurement state is 
diagonal 



t—s = n — h 
t — s = n — k. 



We prove that the post-measurement state 

q£ = Tr s [H ® \n)(n\ g 1 ® \n)(n\}/p n (Al) 

of STS and MTS after local measurement in the number 
basis is diagonal (here, p n = Tr[pl <E) \n)(n\]). We have 



and this implies h = k; therefore the post- measurement 
state is diagonal in the number basis: 

(<£)M = S lhk J2p t s h (N 1 )pt h (N 2 )\O hk ( S t)\ 2 (A7) 
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Appendix B: Discord does not depend on the phase 
of displacement 

We show that the (non-Gaussian) discord in the 
displaced number basis does not depend on the 
phase of displacement for STS and MTS. The ar- 
guments is best given in the characteristic func- 
tion representation of the states citeParisbook The 
STS and MTS states have a Gaussian characteris- 
tic function x[g](A) = exp(— ^A T aA) where A = 
^= (ReA^ , IiyiXa , Re As , Im As ) and the covariance matrix 
is given by 



A C 
C T B 



a 
a 



c 
±c 



c b 
±c 



(Bl) 



where ±c is — c in the case of STS and 
MTS. For STS we have 



-c in the case of 



X[q](\a,\b) = exp(-a|A A | 2 - 6|A B | 2 + 2cRe[X A \ B }) 

(B2) 

while for MTS the same expression holds upon changing 
2cRe[AAAs] — > 2cRe[A^As]. In the following, we shall 
carry on the argument for STS, but the MTS case is 
fully equivalent. If we perform a displacement on one 
mode, g —> D(a)gD^ (a) = g a , the effect on the charac- 
teristic function is easy to evaluate. Using the relation 
D(a)D(X)D^(a) = D(X) exp(— 2iIm[Aa*]) {33J we obtain 

x[Qa](^A, As) = x[£>](Aa, X b ) exp(-2iIm[A_BQ;*]) (B3) 

Suppose we perform a masurement on mode B in the 
number basis {1I„ — \n)(n\}. The post-measurement 
state of mode A is g^ n — ^-TrB[£> Q II n ] where p n = 
Tr[g a Jl n ]. By use of the trace formula [33[ 

TrpiOs] = - [ d 2m X x[Oi}(X) X [0 2 }(-X) 

7T J C m 

we obtain the characteristic function 

x[giJ x A) = — I d 2 x B x [n„](A B ) x [Q a ]{x A ,-\ B ) 

(B4) 

Since x[n n ](As) = e~5l A sl 2 L„(|A B | 2 ), where L n is the 

Laguerre polynomial L n {\X B \ 2 ) = E"=o L~t) l ^ir~ > we 
have explicitly 

x[0oJ(Aa) = — / d 2 X B L n (|A B | 2 )exp(-a|AA| 2 

KPn Jc 

-(b + 1/2)|A B | 2 - 2cRe[AAA B ] - 2iIm[\ B a*]) 



(B5) 

In order to see that this expression depends on |a| only 
we can implement the change a — > a' = e l8 a and we 
have 

x [g^}(X A ) = — [ d 2 X B L„(|A B | 2 )cxp(-a|A A | 2 

KPn Jc 

-(b + 1/2)|A B | 2 - 2cRe[X A X B ] - 2ilm[X B a* e- 10 ]) 



By changing variable X B — > e l8 X B we see that 

X [gi.JX A ) = X [gi n ]^Ae- ld ) (B6) 

Therefore, we have xiQa'J ^ xbaj) hence gfc, ^ g£ n . 

However, g A , and g A n have the same spectrum. Indeed 

x[Qa n (^A)] and x[Q A n (^Ae l6 )] are related by a simple 
quadrature rotation 



(B7) 
(B8) 



a\ — > a\e 
which means that 

gi n =Ug A U^ 

where U is the free evolution of mode A, U = e i6a ' iai . 
Since g A , and g A n are related by a unitary, they have 
the same spectrum. Therefore, the spectrum (hence, the 
entropy) of the reduced post-measurement state g A does 
not depend on the phase of a but just on |a|. If follows 
that the non-Gaussian quantum discord of g a does not 
depend on the phase of a, QED. 

As for the non-Gaussian geometric discord, it is obtained 
as 

Tr[(g) 2 ]-^Tr[gl^®nf] = 

n 

= Tr[(g) 2 ]-]TTr[^Xj 

n 

By the same arguments before, leading to eq. (|B8j) . we 
immediately see that the second trace does not depend on 
the phase of a, hence the geometric discord does neither. 



Appendix C: Understanding the behavior for 
growing a 

Let us now consider in detail the bahaviour for growing 
a. We will show that the non- Gaussian discord in the 
displaced number basis tends to the Gaussian discord 
as the displacement tends to infinity, — > D Q as 
a — > oo. 

First, we will show that 



as a 



oo 



(CI) 



This is best shown in the characteristic function formal- 
ism. The post-measurement state of mode A has the 
characteristic function (|B5|) . Since the phase of a is irrel- 
evant for the discord, we will assume a € K in the follow- 
ing. The post-measurement state characteristic function, 
Eq. (|B5[) . is the Gaussian integral of a polynomial. By 
using a well-known trick of Gaussian integrals, we can 
rewrite 

X^IKAa) = — / d 2 X B L n (d/drj) x 

exp (—7 |A B | 2 + 2cReA j4 ReA B - (2cIm[A, 4 ] + 2ia)Im[A B ] 
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(a) Re 




(b)Im 

FIG. 9: x[Qc,a]W for a = 0.1 (red), a = 1 (green), a = 5 
(blue) and n = (solid), n = 1 (dashed), n = 2 (dotted). 
Here, we have A = (N T + 1/2) cosh A, C = (N T + 1/2) sinh A 
with A = 0.5, N T = 0.5. 



where 7 = 6 + 1/2 and the formal expression L n (d/dy) 
means J^iLo hw*' This expression can now be 
moved outside the integral, so that we are now left with 
a purely Gaussian integral of the form 



d 2 A H e -^ T sMA B+ AlB 



where M = diag{47,47}, B = (2cReAyi, -2cImAyt + 
2ia), Ab = (ReA_B,ImAs). The integral gives 

2x — e \B T M- 1 B so that we fi na ii y get 
Vdetjvi 

x[et\^A) = -e~i a ^ 2 L n (-d/dj) x (C2) 

Pn 

1 ^ ^ /^c 2 |Aa| 2 - a 2 ~ 2iacImA^ 



— exp 

7 



27 



Let us define x = c 2 \Xa\ 2 — a 2 — 2iaclm\A- Then we 
have 

.1 /c 2 |A A | 2 - a 2 - 2iacImA J 4 \ 
L n (-d/d7)-exp I — \ = 

= L n (-d/d 7 )-e x ^ = F n fr,z)e?'*> 
1 

where F n (^ : x) = J^k fk{l) xk is necessarily a polyno- 
mial of degree n in with 7-dependent coefficients fk{l)- 
Therefore, 



xbl](AA) = -e- (Q - c2(b+1/2rl)l ^ |2 x 

Pn 

e -zc( 6+ l/2)-WmA Ae - a V27 Fn(xx ) 



(C3) 



The norm is 



= e -(a-c 2 (fc+l/2)- 1 )|A A | 2 e - lC (fc+l/2)- 1 a/ m A A x 



- a2/27 Fn(7^) 



Aa=0 



- 2 /27 Fn( _ a2) 



so that 



-ic(b+l/2)- 1 aIm\ A Frijl, x) 

^„(7,-« 2 ) 



(C4) 



- _™2 



This function is exponentially decaying as e~ s ^ XA ^ where 
s = a — c 2 / (6+1/2), hence it is vanishing for | A^4 1 2 3> 1/s. 
Therefore, we can consider values of \\a\ 2 in the region 
\ 2 A < 1/s. In this region, we we have lim Q _j 
because a^$> Xa and thus 

li F n(l,x) _ fnh)a 2 n = 
a^oo F„(7, -a 2 ) f n (l)a 2 n 

In conclusion, as a — > 00 we have 



(C5) 



which implies the desired result (jCip . Q-E-D. 

This result means that the conditional state of A is 
independent of n and equal to the n — result. In fig. 
EH we show x[fi ) Q ],x[£'Q 1 ]:X[£ , a 2 ] for growing values of a. 
The three curves converge already for a ~ 5. As a con- 
sequence of n-independence, we have 



?n,A/"g 



(4B) = yVsG&)->sG4) 



(C6) 



But p„ is just the post-measurement state correspond- 
ing to POVM element D(a)\0){0\D^ (a) = \a)(a\, i.e, 
a Gaussian state with covariance matrix up = A — 
C(B + I/2)~ 1 C T (Schur complement), and mean fip = 
X(B + l/2)- 1 C T , where X = (a,0). On the other 
hand, from the discussion in sec. HI] we know that the 
optimal Gaussian POVM is a heterodyne measurement 
{IL3 = D(fi)\Q){Q\D*{fi) = \P){P\}. In this case, as al- 
ready explained in sec. HU the entropy of the post mea- 
surement state g"p is independent of the measurement 
result /3, hence the conditional entropy coincides with 
the entropy of of the (3 — a result. Therefore, we also 
have S g (A\B) = S(g Q0 ). Therefore, we conclude that 
the non Gaussian discord D Mg {A : B) in the displaced 
number basis tends to the Gaussian discord D g (A : B) 
as a — > 00, QED. 

To be rigorous, we did not prove that the D Mg {A : B) 
is lower bounded by D g (A : B), and we cannot rule 
out the possibility that D Mg (A : B) < D g {A : B) 
for intermediate values of a. However, our numerical 
data do not support this possibility since we never ob- 
serve D Ng {A : B) < D G (A : B) and we expect that 
D Ng {A : B) -> D g (A : B) fr om above as a — > 00. 
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